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We report measurements of the magnetic excitation spectrum of the layered antiferromagnet 
La2Co04 by time-of- flight neutron inelastic scattering. In the energy range probed in our experi- 
ments (0-250 meV) the magnetic spectrum consists of spin-wave modes with strong in-plane disper- 
sion extending up to 60meV, and a nearly dispersionless peak at 190 meV. The spin-wave modes 
exhibit a small (~ 1 meV) dispersion along the magnetic zone boundary. We show that the magnetic 
spectrum can be described very well by a model of a Heisenberg antiferromagnet that includes the 
full spin and orbital degrees of freedom of Co 2+ in an axially-distorted crystal field. The collective 
magnetic dynamics are found to be controlled by dominant nearest-neighbour exchange interactions, 
strong XY-like single-ion anisotropy and a substantial unquenched orbital angular momentum. 

PACS numbers: 75.50.Ec, 75.30.Et, 75.30.Ds, 78.70.Nx 



I. INTRODUCTION 

The evolution from antifcrromagnctism to high- 
temperature superconductivity with carrier doping of the 
layered copper oxides^ has inspired a vast literature on 
the electronic, structural, dynamical, and chemical prop- 
erties of related materials. It has become clear from 
these studies that superconductivity and commensurate 
antifcrromagnetic order arc only two out of many differ- 
ent competing ordering tendencies found in systems of 
strongly interacting electrons. 

Among other forms of order found in doped Mott in- 
sulators are nematic phases characterised by unidirec- 
tional density- wave states involving combined charge and 
spin order. Such 'striped' phases were first identified 
many years ago in hole-doped (La,Nd)2CuC>4 (Ref. 0) 
and La2NiC>4 (Refs.0-@), but their significance for high- 
temperature superconductivity has been the subject of a 
continuing debate. Although much of the focus has been 
on the cuprates, the nickelates have contributed to this 
debate on account of their relatively well correlated and 
stable stripe order which is amenable to experimental 
investigation. One drawback, however, is that holes lo- 
calised on Ni 3+ ions in hole-doped L^NiCU carry a spin 
which can interact magnetically both with other spins 
in the charge stripes and with the surrounding antifcr- 
romagnetic matrix of Ni 2+ . The influence of these in- 
teracting magnetic degrees of freedom on the properties 
of stripes in nickelates has yet to be fully evaluated, but 
spin correlations associated with both Ni sites have been 
observed^ and there remain some unexplained features in 
the spin excitation spectral 

Recently, evidence has been presented for the exis- 
tence of stripe phases in the layered cobaltate system 
Laa-xSr^CoCU (Ref.|), which is isostructural with hole- 
doped La2Cu04 and La2Ni04. Neutron diffraction mea- 
surements on half-doped (x = 0.5) L^CoCU show clear 
evidence for a charge ordering of Co 2+ and Co 3+ ions in 
a checkerboard pattern below T co « 825 K, with an an- 



tifcrromagnetic ordering transition at a much lower tem- 
perature Tn ~ 60Kii£~— The antiferromagnetic order 
in La2- K Sr z: Co04 in the range 0.3 < x < 0.6 was ob- 
served to be periodically modulated with the modulation 
wavevector nearly proportional to a characteristic of 
the stripe phase found in La2- x Sr 2 .Ni04. The Sr-dopcd 
cobaltate system has an advantage over the nickelates in 
that for the compositions in which stripe-like magnetic 
correlations occur there is evidence to suggest that the 
Co 3+ ions adopt the low spin (S = 0) state at low tem- 
peratures and are therefore not magnetically active*^— 
Hence, the layered cobaltates offer the chance to inves- 
tigate the fundamental interactions and excited states of 
an ordered stripe phase in which the doped holes do not 
possess low-energy spin degrees of freedom. 

Attempts to understand the electronic phases in 
L^-zSr^CoCU will require some basic knowledge of the 
parent antiferromagnet La2Co04. Although the crys- 
tal structure and magnetic order of La2Co04 have been 
studied in detail^ no measurements of the magnetic ex- 
citation spectrum have been published before now. 

In this paper we present high-resolution measurements 
of the magnetic spectrum of L^CoC^+s (<5 w 0) by neu- 
tron inelastic scattering. The methodology relates closely 
to that employed in a recent investigation carried out on 
the half-doped cobaltate Lai.sSro.sCoC^ (Ref. [TH). The 
data presented here on La2Co04 extend up to 250 meV 
in energy and reveal spin-wave-like excitations with a 
bandwidth of 60 meV. At much higher energies, there 
are excitations of a localised character. A good descrip- 
tion of the data is achieved with a spin-wave model for a 
quasi-two-dimensional antiferromagnet that includes the 
full spin and orbital degrees of freedom of the Co 2+ ions. 
The results show that La2Co04 has dominant nearest- 
neighbour exchange interactions, although a weak disper- 
sion along the zone boundary indicates that more distant 
interactions or non-linear terms in the Hamiltonian are 
not negligible. The anisotropy is strongly XY-like but 
there is also a weak in-plane anisotropy. 
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FIG. 1. (a) In-plane magnetic structure of La2Co04. The 
dashed square shows the conventional /4/mmm unit cell of 
the HTT phase, and the filled square represents the magnetic 
unit cell, which coincides with the y/2 x y/2 chemical unit cell 
of the LTT phase. The exchange interactions used to model 
the magnetic spectrum are indicated, (b) Diagram of the re- 
ciprocal space lattice corresponding to the I4/mmm cell. The 
filled square indicates the magnetic Brillouin zone centered on 
(0.5,0.5). The dashed lines show the path through recipro- 
cal space along high-symmetry directions used for detailed 
analysis of the magnetic excitation spectrum — see Fig. [5] 



II. CRYSTAL AND MAGNETIC STRUCTURE 

In common with the stoichiometric La2Cu04 and 
La2NiC>4 compounds, La2CoC>4 exhibits three different 
structural phasesJ^ (i) T > T\ high-temperature tetrag- 
onal (HTT), space group IA/mmm] (ii) T 2 < T < 
T\ low-temperature orthorhombic (LTO), space group 
Cmca; (iii) T < T 2 low-temperature tetragonal (LTT), 
space group PA 2 /ncm. The structural transition tem- 
peratures for La 2 Co04 are T\ « 900 K (Ref. 0) and 
T 2 = 120 - 135 K (Ref. [H). The latter is reported to 
be first order. Throughout this paper we shall use the 
conventional IA/mmm unit cell as a basis for the re- 
ciprocal lattice. The low temperature lattice constants 
referred to this cell are a = b = 3.91 A and c = 12.6 A. 
The true LTT unit cell has in-plane dimensions which are 
\[2 x \/2 larger than those of the IA/mmm pseudo-cell 
— see Fig. [TJa). 

The transition to magnetic order occurs at Tn ~ 275 K, 
and a magnetic reorientation occurs at T 2 coincident 
with the LTO-LTT structural transition. The antifcr- 
romagnctic structure has an ordering wavevector q m = 
(0.5,0.5,0), with ordered moments lying in the C0O2 
plane. Assuming collinear order, the difference between 
the magnetic structures in the LTT and LTO phases is 
that in the LTT phase the moments are perpendicular 
to q m whereas in the LTO phase they are parallel to 
q m . The distinction between these structures depends 
on the relationship between adjacent layers. Another 
possibility is that the structure is collinear within the 
layers but the moment direction rotates by ±90° from 
one layer to the nextJ£ In the absence of inter-layer cou- 
pling all these structures have the same energy. As we 
did not observe any evidence in the excitation spectrum 
for inter-layer coupling we will treat the magnetic order 
as two-dimensional. Figure da) shows the in-plane mag- 



netic order with the moments arbitrarily chosen to point 
along the horizontal axis. 



III. EXPERIMENTAL DETAILS 

A single- crystal sample of La 2 Co04 with a mass of 
approximately 5 g was grown in Oxford by the optical 
floating-zone method. Polycrystalline La 2 Co04 was pre- 
pared from high-purity (>99.99%) La 2 03 and C03O4 by 
solid-state reaction. Stoichiometric amounts of the oxides 
were mixed and reacted at 1050°C for 48 hours under a 
flowing atmosphere of CO/C0 2 mixed in the ratio 1:10. 
A reducing atmosphere is needed to avoid the formation 
of LaCo03. The powder was re-ground and sintered at 
1100°C in a flow of argon for a further 48 hours. No im- 
purity phases could be detected in the product by x-ray 
powder diffraction. The La 2 Co04 powder was isostat- 
ically pressed into rods of diameter 12 mm and length 
120 mm. The rods were sintered in an argon atmosphere 
at 1250°C for 24 hours. Crystal growth was carried out 
in a four-mirror image furnace (Crystal Systems Corpo- 
ration) in flowing argon at a growth speed of approxi- 
mately 2mmhr~ 1 with counter- rotation of the feed and 
seed rods at 25 rpm. 

Crystals grown by this method contain an excess of 
oxygen. To achieve stoichiometry the as-grown crystal 
was annealed at 850°C for 72 hours in flowing CO/CO2 
(1:10 ratio). A fragment of the annealed crystal was 
ground to a powder and subjected to a thermogravimet- 
ric analysis. From the measured weight loss we deter- 
mined the oxygen nonstoichiomctry to be 6 = —0.03 ± 
0.02. This suggests that the crystal is close to the ideal 
stoichiometry, if anything slightly oxygen-deficient. 

Magnetisation measurements were performed with a 
superconducting quantum interference device (SQUID) 
magnetometer (Quantum Design) on a small crystal cut 
from the same rod as the neutron scattering crystal. 
Measurements were made by the dc method with a mea- 
suring field of strength 1000 Oe (fioH = 0.1 T) applied 
along the crystallographic c axis. Zero-field-cooled (ZFC) 
data were recorded on warming after the sample had been 
initially cooled from 350 K in zero applied field, and field- 
cooled (FC) data were recorded while cooling the sample 
from 350 K in the measuring field. 

Unpolarised-ncutron inelastic scattering measure- 
ments were performed on the direct-geometry chopper 
spectrometer MAPS at the ISIS facility^ Neutron time- 
of-flight instruments with large position sensitive detec- 
tor arrays such as MAPS allow sampling of vast regions of 
(Q, lu) space simultaneously, where Q and fiw are, respec- 
tively, the wavevector and energy transferred from the 
neutron to the sample. This is very advantageous in stud- 
ies where the excitation spectrum is required throughout 
the Brillouin zone. 

In preparation for the neutron measurements the 
La2Co04 crystal was sealed in a thin- walled aluminium 
can containing helium exchange gas and aligned with the 
c axis parallel to the direction of the incident-neutron 
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beam. Cooling was provided by a closed-cycle refrigera- 
tor. Data were collected with incident-neutron energies 
of 51, 86, 111, 152 and 303 meV. The energy resolution 
was typically 5% of the incident energy (full width at half 
maximum) at zero energy transfer, decreasing slightly 
with increasing energy transfer. Under the chosen experi- 
mental conditions the wavevector resolution is largely de- 
termined by the divergence of the incident-neutron beam 
which is approximately 0.5°. Spectra from La2Co04 were 
recorded at several temperatures between 6 and 300 K. 
Separate measurements of a standard vanadium sample 
were made at each incident energy to normalise the spec- 
tra and place them on an absolute intensity scale. 

For presentation and analysis, the neutron data were 
transformed from raw time-of-flight spectra into an in- 
tensity map as a function of Q and hu>. With a fixed 
sample orientation, only three out of the four compo- 
nents of (Q,w) are independent. We chose the two in- 
plane wavevector components (Q x ,Qy) = (h,k) x 2n/a 
and energy as the independent variables, which means 
that the out-of-plane wavevector Q z = I x 2ir/c varies 
implicitly with energy transfer. For a two-dimensional 
scattering system, however, there is no dispersion in the 
out-of-plane direction and the gradual variation of scat- 
tering intensity with Q z can be included in a model (and 
was done so in this work). The justification for treat- 
ing La2CoC>4 as a two-dimensional magnetic system is 
that the magnetic spectra show no discernible periodic 
modulation in intensity with Q z (i.e., with fkS). 

In order to quantify the magnetic dispersion we 
performed a series of constant-energy and constant- 
wavevector cuts through the data volume along high- 
symmetry directions using the MSLICE software^ Be- 
fore performing these cuts, data at symmetry-equivalent 
wavevectors were averaged to improve the signal. 



IV. RESULTS 

The temperature dependence of the FC and ZFC sus- 
ceptibility (x = M/H) is shown in Fig. [2] (inset). Both 
curves show a change in slope at approximately 276 K 
consistent with the antifcrromagnctic transition, and 
sharp anomalies at 124 K close to the temperature T2 
at which the LTO-LTT structural transition is expected. 
The onset of antiferromagnetism at Tn ~ 276 K is con- 
firmed by the temperature dependence of the neutron 
diffraction intensity recorded at the magnetic Bragg peak 
position (1.5, 1.5, 1), shown in Fig. [2j To estimate Tn we 
fitted the data to a power law J oc (1 — T/Tn) 2 ^ , assum- 
ing a Gaussian distribution of Neel temperatures about 
the mean value (Tn) with standard deviation of or- This 
function was found to give a good description of diffrac- 
tion data near Tn in Ref. [l6l . The parameters obtained 
from our data were (Tn) = 275.5(5) K, <jt = 1(1) K, and 
j3 = 0.15(1). The transition temperatures measured on 
our sample are consistent with previously reported values 
of Tn = 275 K and T2 ss 135 K for a nominally stoichio- 
metric crystal^ Although there is some discussion in 
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FIG. 2. Main figure: temperature dependence of the 
(1.5, 1.5, 1) magnetic Bragg peak of La2Co04 showing the on- 
set of magnetic order at Tn ~ 276 K. The solid line shows a 
power law fit to data with a Gaussian distribution of Neel 
temperatures. Inset: field-cooled (FC) and zero- field-cooled 
(ZFC) magnetic susceptibility of La2Co04 recorded with a 
measuring field of 1000 Oe applied parallel to the c axis. The 
magnetic ordering transition and the LTO-LTT structural 
transition are indicated by arrows. 



the literature about the precise composition of La2CoC>4 
prepared under different conditions ; 16 > 19 ~ — we can at the 
very least be confident that our crystal is close in com- 
position to the one used in Ref. [TH 

We note that the FC and ZFC susceptibility curves 
separate below 350 K, which is not expected in the para- 
magnetic phase. This indicates that the sample contains 
a small amount of ferromagnetic impurity. The FC-ZFC 
separation was not observed in the as-grown crystal. The 
most probable explanation is that a tiny amount of el- 
emental Co was formed during the CO/CO2 annealing 
step. This is consistent with the slight oxygen deficiency 
found from the thermogravimetric analysis. As there is 
no unexplained secondary signal in our neutron scatter- 
ing spectra this impurity must be present in very small 
quantities so is of no consequence to our neutron results, 
but it does mean that the susceptibility curves shown in 
Fig. [2] contain a background signal in additional to the 
signal from pure La2Co04. 

We now turn to the neutron scattering spectra. Fig- 
ure [3] provides an overview of the data collected at 6K. 
Panels (a)-(c) are constant-energy slices at three differ- 
ent energies, and panels (g)-(i) are energy-Q slices to 
illustrate the magnetic dispersion. The spectrum is dom- 
inated by a spin-wave-like conical dispersion which rises 
from the in-plane antifcrromagnctic ordering wavevector 
q m = (0.5, 0.5) and equivalent positions [the M-points of 
the square-lattice Brillouin zone — see Fig. UF>)]. This 
mode has a gap of approximately 10 meV at the M-point 
and rises to a maximum energy of 60 meV at the E-point 
on the Brillouin zone boundary. A much weaker branch, 
displaying an upwards dispersion with a minimum energy 
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(a) E = 10-12 meV (d) E = 10-12 meV 
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FIG. 3. Measured and simulated magnetic spectra of La2CoC>4 at 6K. Panels (a)-(c) show intensity maps averaged over 
2meV energy ranges centered on different energies as indicated, with the corresponding calculated spectra shown in (d)-(f). 
The magnetic dispersion along two high-symmetry directions is displayed in (g)-(i), with corresponding simulations in (j)-(l). 
Data in (a)-(h) were measured with an incident neutron energy Ei = 86meV, while (i) was measured with Ei = 303 meV. 
The units of intensity indicated by the colourbars are mb sr _1 meV -1 f.u. -1 . The simulated spectra are calculated from the 
spin-orbital spin- wave model described in the text. 
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FIG. 4. Constant-Q cuts taken through the measured and 
simulated spectra. Panel (a) shows a series of cuts follow- 
ing the dispersion from the zone center (0.5, 0.5) to the an- 
tiferromagnetic zone boundary at (0.75,0.75) measured with 
incident energy lllmeV. Panel (b) shows the spin gap of 
~10meV at (0.5,0.5) measured with an incident energy of 
51 meV which gave an improved energy resolution. Panel 
(c) shows the measured and simulated excitation mode at 
~ 190 meV. These data were obtained with Ei — 303 meV. 



at M of 46 meV, corresponds to the first mode translated 
by q m . The large splitting of the modes at M shows that 
the anisotropy is strongly XY-like. The lower and upper 



modes correspond to in-plane and out-of-plane fluctua- 
tions, respectively. Figurc|30) shows data up to the max- 
imum energy explored in our experiment. This reveals 
only one other significant feature — a band of scatter- 
ing in a narrow range of energies close to 190 me V. We 
note that scattering from phonons is much weaker than 
magnetic scattering in the range of Q studied. 

To give a different impression of the data we present 
in Fig. |4] examples of constant-Q cuts taken through 
the data volumes measured with incident energies Ei — 
51, 111 and 303 meV. To extract the magnetic disper- 
sion in a form suitable for fitting to a model we per- 
formed a large number of such constant-energy cuts at 
wavevectors along the reciprocal-space paths indicated 
in Fig. []Jb). The peaks in these plus some additional 
constant- wavevector cuts were fitted with Gaussian func- 
tions on a linear background. The peak centers deter- 
mined this way are plotted along the high symmetry di- 
rections in Fig. [S] 

An interesting behaviour is observed along the mag- 
netic zone boundary: the energy of the magnon branch 
is not constant but varies by approximately 1.5 meV. As 
discussed below, this is significant because a dispersion 
along the zone boundary indicates a need to go beyond 
a linear spin-wave model with nearest-neighbour interac- 
tions only. To emphasise this effect, we show in Fig.^a) 
the energy and (b) the integrated intensity of the magnon 
peak along the entire length of a zone boundary (XSX). 
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FIG. 5. The lower figure shows the dispersion of the mag- 
netic excitations of La2CoC>4 along high symmetry directions 
in the 2D Brillouin zone defined in Fig. [TJ Open circles are 
points extracted from cuts through the measured data vol- 
ume. The lines show the dispersion of the modes calculated 
with the many-level spin-wave model described in the text. 
The upper figure shows the response functions S ac " for each 
mode calculated from the many-level model. The normaliza- 
tion of the response functions is per formula unit of La2Co04. 



The maximum in the dispersion at £ is seen to coincide 
with a minimum in its intensity. Because the disper- 
sion surface forms a ridge along the zone boundary care 
was taken to select an appropriately-sized box in Q over 
which to average the data so as to avoid systematic er- 
rors from the curvature of the dispersion surface while at 
the same time having good enough statistics to extract 
the peak energies and integrated intensities. Figure 
shows energy cuts taken at an X-point and a E-point to 
illustrate the difference between the magnon peaks at the 
zone corner and zone edge. 



V. ANALYSIS AND DISCUSSION 

Magnetism in Co 2+ compounds such as La2Co04 
is generally influenced to a significant degree by un- 
quenched orbital angular momentum which is respon- 
sible for, among other things, the strong anisotropy in 



the susceptibility observed in many such compounds. In 
a recent studyi^ of the magnetic excitations in the half- 
doped cobaltate Lai.sSro.sCoG^, which is also an antifer- 
romagnet, we developed a model to describe the magnetic 
spectrum including both the spin and orbital angular mo- 
mentum of the Co 2+ in the high-spin configuration (3d 7 , 
S = 3/2, L = 3). The model is an advance over con- 
ventional (spin-only) spin- wave theory in that it includes 
level-mixing within the 2 +1 L term caused by the ligand 
and exchange fields, and hence the parameters that de- 
scribe the single-ion anisotropy and exchange interactions 
are physically realistic. As far as the magnetic spectrum 
is concerned, the admixture of basis states means that ex- 
citations to levels above the first excited single-ion level 
can propagate and can be observed by neutron scatter- 
ing. Moreover, the orbital component of the single-ion 
states needs to be included for an accurate calculation of 
the neutron cross section. 

The model employs the Hamiltonian 

H = 52J ik S r S k 



E 



^/;nr:L ; < • ai.,.s ; • n; s 



l.m 



•(i) 



The first term describes an isotropic Heisenberg exchange 
interaction between pairs of S = 3/2 spins. For La2Co04 
we include only the nearest-neighbour and next-nearest- 
neighbour exchange interactions J, J\ and J2, as defined 
in Fig.[TJ The remaining terms in ([1} arc single-ion terms. 
The first of these represents the crystal (ligand) field act- 
ing on the Co 2+ ions. The Of 1 are Stevens operator- 
equivalents with B™ the corresponding crystal-field pa- 
rameters. The axially-distorted octahedral crystal field 
from the neighbouring 2 ~ ions is described by the op- 
erators O2, O4 and Of. We kept the same values for 
the parameters B® and B\ as found for Lai.sSro.sCoO^ 
in Ref. OJ: B\ = -1.35 mcV and B% = -8.00 meV. 
These are estimated from a point-charge calculation and 
scaled to match the cubic crystal field splitting observed 
in CoO.— The parameter B\ controls the out-of-planc 
anisotropy and was adjusted to obtain a good fit to the 
magnetic spectrum. Its final value (see below) differs 
from that deduced for Lai. 5 Sr .5CoO4 by only ~10%. 
The term AL • S is the spin-orbit coupling. The cou- 
pling constant A = — 18.7meV used here has been de- 
duced from reflectivity measurements of CoO by optical 
spectroscopy^ The final term H a • S represents a small 
uniaxial anisotropy which defines the in-plane orientation 
of the moments and produces a spin gap at the T-point 
(and, equivalently, the M-point). We chose the moments 
to lie along the x axis, and to achieve this the anisotropy 
field H a points along +x on one of the antifcrromagnctic 
sublatticcs and along —a; on the other. 

The partial differential scattering cross-section de- 
pends on the response functions ^""(Q, u>) describing act 
magnetic correlations. In the dipole approximation the 
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relation is^ 

a 

(2) 

where 

S aa (Q,u) = IOW a (Q)|0)| 2 5[w - Wj(Q)]- (3) 

3 

Here, fc; and kf are initial and final neutron wavevectors, 
(jr /2) 2 = 72.8 mb, f(Q) is the dipole magnetic form 
factor of Co 2+ , e~ 2W is the Debye- Waller factor which 
is close to unity at low temperatures, and Q a = Q Q /|Q| 
is the a component of a unit vector in the direction of 
Q. The response function (per La2Co04 f.u.) described 
in Eq. ([3]) takes into account both the spin and orbital 
magnetization M = — (L + 2S) in the transition matrix 
element connecting the ground state to an excited mode 
j. The procedure to diagonalise the Hamiltonian (JTJ to 
obtain the dispersion and response functions of the mag- 
netic modes is described in detail in Ref. [lH 

The parameters of the model were refined from a fit 
to the measured dispersion carried out by a simulated- 
annealing algorithm. Because of the computer time re- 
quired to diagonalise the Hamiltonian for the complete 
set of 2 x {(2L + 1)(25 + 1) - 1} = 54 excited states 
(twice the number of single-ion excited states because we 
have two magnetic sublattices) we restricted the num- 
ber of observables included in the fit to just enough to 
represent all the important features of the data, includ- 
ing the high-energy signal at ~190mcV. The parame- 
ters varied in the fit were -6°, J, Ji, J-2, and H a . The 
best fit was achieved with parameters B\ = 14.6(1) meV, 
J = 9.69(2) meV, J x = 0.14(2) meV, J 2 = 0.43(1) meV, 
and fP = 0.66(6) meV. The calculated dispersion and 
response functions of the magnetic modes are shown in 
Fig. [5] together with the full set of data points for the 
lowest energy modes determined from the measurements. 
The agreement is seen to be very good. The fit indicates 
that the next-nearest-neighbour exchange constants J\ 
and J2 are very small but not zero. As a test, we re- 
peated the fit with Ji and Ji fixed to zero and found that 
the quality of best fit worsened, as indicated by the stan- 
dard goodness-of-fit parameter x 2 per degree of freedom 
which increased from 4.5 to 11.1. Therefore, we believe 
that the obtained values of J\ and J2, though small, are 
significant. 

To further visualise and assess the model we calcu- 
lated intensity maps and cuts to simulate those obtained 
from the experiment. Figures [3] and [4] show the simu- 
lations alongside the corresponding experimental data. 
The quantity plotted is (ki/kf)d 2 a/d£ldE{ per formula 
unit (f.u.), i.e., the partial differential cross section mul- 
tiplied by a factor k\jk{ as defined in Eq. ([2])- The dipole 
magnetic form factor of Co 2+ and the Q orientation fac- 
tor that determines the weighting of the different re- 
sponse functions are included in the simulated spectra. 
The simulations also take into account a number of other 
experimental factors: (i) we averaged over a 50:50 mix- 
ture of equivalent magnetic domains in which the ordered 
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FIG. 6. Dispersion of the magnon peak along the magnetic 
zone boundary in La2CoC>4. Variation of (a) the peak po- 
sition and (b) the integrated intensity of the magnon peak 
in constant-Q cuts, (c) Constant-Q cuts at Q = (0.75, 0.75) 
and (0.5, 1) fitted with Gaussian lineshapes. The data are 
from the run with incident energy lllmeV and sample tem- 
perature 6 K. 



moments point along the x and y axes, respectively; (ii) 
the spectra are broadened in energy and wavevector by 
the estimated resolution of the MAPS spectrometer (see 
Sec. IIII[) : (iii) we included an estimate of the absorption 
and self-shielding of the neutron beam by the sample, 
which reduces the intensity by a factor of typically 0.65- 
0.80 depending on the incident neutron energy and fvjj. 
An additional scale factor of 0.4 was applied uniformly 
to all calculated spectra in order to match the measured 
absolute scattering intensity. 

The simulations show that the model provides a very 
good description of the entire observed spectrum of 
La2Co04. The relative intensities of the magnetic exci- 
tations are reproduced to within 10-20%, including the 
band of scattering at ~190meV, which from Fig.[5]is seen 
to originate from a mode with longitudinal (xx) char- 
acter together with some less-intense transverse modes. 
Magnetic excitations are also present in the model at 
around ~115 and ~165meV but are predicted to carry 
negligible spectral weight and are not observed — see 
Fig. IU[c). The additional scale factor of 0.4 needed to 
match the absolute intensity is similar to that required 
for Lai.5Sr .5CoO4 (Ref. [l5|). It is accounted for partly 
by the size of the ordered moment. The observed ordered 
moment is 2.9 /xb (Ref. Htl ) whereas the ordered moment 
in the (ionic) model is 4.1 /ib- The difference between 
observed and calculated moments may be an effect of 
covalency, which would also modify the magnetic form 
factor relative to the free ion form factor in such a way 
that could cause an additional reduction in intensity, as 
recently found in a cuprate chain compound* 2 ^ 

For reference, we also compared the low-energy part 
of the spectrum (hui < 60mcV) with standard linear 
spin-wave theory for an effective spin-i antiferromag- 
net, which neglects the orbital component of the modes. 
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We used the same model as described in Ref. [jja in 
which the magnetic anisotropy is described by anisotropic 
nearest-neighbour exchange interactions J x = J(l + e), 
J y = J, and J z = J(l — 5). The parameters e and 8 
control the in-planc and out-of-plane anisotropy, respec- 
tively. The more distant interactions J\ and Ji were 
included too, but because they are relatively small we 
treated these as isotropic. Wc found that the lower en- 
ergy modes can be well described by this model. In fact, 
an equally good description of the data (as reflected in 
the value of % 2 ) could be found with sets of parameters 
in which J\ and Ji are both positive or both negative: (i) 
J = 9.89(1) meV, Ji = 0.04(1) mcV, J 2 = 0.13(1) meV, 
e = 0.013(1), 5 = 0.283(4), or (ii) J = 8.30(6) mcV, 
Ji = -0.35(2) meV, J 2 = -0.63(3) meV, e = 0.024(1), 
S = 0.383(5). By contrast, the spin-orbital many-level 
model clearly favors the case with J\ and Ji both pos- 
itive. The spin-orbital model can discriminate the two 
cases because of the inclusion of the higher excited levels. 
Only the parameter set with J\ and Ji both positive fits 
the low energy modes (E < 60 meV) and reproduces the 
peak in the spectrum at ~190 meV and absence of any 
other measurable peaks between 60 and 250 meV. An- 
other drawback of the effective spin-^ linear spin- wave 
model is that the intensities are not accurately described 
because of the neglect of the orbital degrees of freedom. 

It is interesting to compare the magnetic spectrum 
of La2Co04 with that of other two-dimensional, square- 
lattice, antifcrromagnctic insulators, particularly in rela- 
tion to the anomalous dispersion along the zone bound- 
ary. By "anomalous" , we mean that the zone-boundary 
dispersion cannot be described within the framework 
of an antiferromagnetic spin-wave model in the linear 
approximation with only nearest-neighbour interactions. 
Inclusion of (i) interactions with more distant neigh- 
bours, or (ii) terms beyond the linear approximation, 
are two ways in which a zone-boundary dispersion can 
be obtained. Other layered antiferromagnets which ex- 
hibit zone-boundary dispersion include La2Cu04 (Refs. 
[H and Hi), Sr 2 Cu30 4 Cl2 (Ref. M) and Cu(DCOO) 2 • 
4D 2 (CFTD, Refs.[28|and[2g|). These are all highly two- 
dimensional, S = I Heisenberg antiferromagnets with 
almost isotropic interactions, and it is thought that the 
zone-boundary dispersion is caused by non-linear terms 
in the nearest-neighbour Heisenberg model. For example, 
in La2Cu04 a model with a four-spin ring exchange was 
employed^ and for CFTD a resonating-valence-bond 
model describing entangled spin-dimer states was pro- 
posed to explain the data^ Interestingly, the behaviour 
along the zone boundary is different in these two materi- 
als: in La2Cu04 both the energy and intensity arc higher 
at X than at E, whereas in CFTD both the energy and 
intensity are higher at E than at X. In La2Co04, on the 
other hand, the energy is a maximum at E while the in- 
tensity is a maximum at X (see Fig. [5]). By contrast, there 
is virtually no zone-boundary dispersion at all in S = 5/2 
square-lattice system Rb2MnF4, Ref. [3tJ In our analy- 
sis of La2Co04 we found that although the zone bound- 
ary dispersion can be satisfactorily reproduced with an 
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FIG. 7. Temperature dependence of the magnetic spec- 
trum of La2CoC>4. The upper panels (a)-(c) show inten- 
sity maps measured at 6, 150 and 300 K along the (h, 0.5) 
direction, (d) Magnetic spectrum at the antiferromagnetic 
ordering wavevector (0.5, 0.5) (M-point) measured at 6, 150 
and 300 K. Inset: temperature evolution of the higher-energy 
magnon mode. 



appropriate choice of J\ and J2, the corresponding in- 
tensity does not have the deep minimum at E found in 
the experiment [Fig. EJb)] . Therefore, whether the zone 
boundary dispersion of La2Co04 is due to interactions 
with more distant spins or arises from quantum effects in 
a non-linear nearest-neighbour model remains an open 
question. 

Finally, we consider the temperature dependence of the 
magnetic spectrum. Figures Ela)-(c) show maps of the 
magnetic scattering measured at T = 6, 150 and 300 K, 
and Fig. [TJd) displays constant-Q cuts at the magnetic 
zone center for the same temperatures. On increasing the 
temperature from 6 to 150 K the 11 meV peak increases in 
intensity due to the increasing thermal population but re- 
mains at the same energy, while the 46 meV peak broad- 
ens and shifts to lower energy [Fig. [TJd) inset] . Although 
La2Co04 undergoes a first-order phase transition coin- 
cident with a magnetic reorientation at T2 w 125 K, the 
in-plane lattice parameters in the LTO phase differ only 
slightly from those in the LTT phase, and the change in 
the magnetic structure only affects the stacking along the 
c axis. It is not surprising, therefore, that the transition 
does not significantly affect the magnetic spectrum. At 
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T = 300 K, the spectrum has become quasielastic and 
there are no longer any sharp inelastic peaks. This indi- 
cates the absence of long-range magnetic correlations for 
T>T N . 



VI. CONCLUSION 

We have measured the excitation spectrum of single- 
crystal La2Co04, an excellent realization of a two- 
dimensional Heisenberg antiferromagnet. We have com- 
bined the experimental results with numerical simula- 
tions to achieve a very good description of the magnetic 
spectrum throughout the entire Brillouin zone, up to an 
energy of 250 meV. The magnetic anisotropy is strongly 
XY-like, but a small uniaxial anisotropy is present which 



will make the low temperature magnetic properties Ising- 
like. An anomalous dispersion along the antiferromag- 
netic zone boundary is observed and can be reproduced 
by including exchange interactions beyond the nearest- 
neighbours but which could also be a manifestation of 
quantum fluctuations in a nearest-neighbour model. 

ACKNOWLEDGMENTS 

We wish to acknowledge helpful discussions with 
B. Keimer, B. Roessli, G. L. Pascut and G. Johnstone 
in the preparation of this work. We thank R. A. Ew- 
ings for valuable assistance with the data analysis, and 
T. G. Pcrring for help with the experiments. P.B. is 
grateful for the provision of a studentship from the UK 
Engineering and Physical Sciences Research Council. 



* peter.babkevich@physics.ox.ac.uk 

1 J. G. Bednorz and K. A. Miiller, Z. Phys. B 64. 189 (1986). 

2 J. M. Tranquada, B. J. Sternleib, J. D. Axe, Y. Nakamura, 
and S. Uchida, Nature (London) 375, 561 (1995). 

3 C. H. Chen, S. W. Cheong, and A. S. Cooper, Phys. Rev. 
Lett. 71, 2461 (1993). 

4 J. M. Tranquada, D. J. Buttrey, V. Sachan, and J. E. 
Lorenzo, Phys. Rev. Lett. 73, 1003 (1994). 

5 P. Wochner, J. M. Tranquada, D. J. Buttrey, and 
V. Sachan, Phys. Rev. B 57, 1066 (1998). 

6 H. Yoshizawa, T. Kakeshita, R. Kajimoto, T. Tanabe, 
T. Katsufuji, and Y. Tokura, Phys. Rev. B 61, R854 
(2000). 

7 A. T. Boothroyd, P. G. Freeman, D. Prabhakaran, 
A. Hiess, M. Enderle, J. Kulda, and F. Altorfer, Phys. 
Rev. Lett. 91, 257201 (2003). 

8 A. T. Boothroyd, D. Prabhakaran, P. G. Freeman, S. J. S. 
Lister, M. Enderle, A. Hiess, and J. Kulda, Phys. Rev. B 
67, 100407(R) (2003). 

9 M. Cwik, M. Benomar, T. Finger, Y. Sidis, D. SenfT, 
M. Reuther, T. Lorenz, and M. Braden, Phys. Rev. Lett. 
102, 057201 (2009). 

10 I. A. Zaliznyak, J. P. Hill, J. M. Tranquada, R. Erwin, and 
Y. Moritomo, Phys. Rev. Lett. 85, 4353 (2000). 

11 A. T. Savici, I. A. Zaliznyak, G. D. Gu, and R. Erwin, 
Phys. Rev. B 75, 184443 (2007). 

12 K. Horigane, H. Hiraka, T. Uchida, K. Yamada, and 
J. Akimitsu, J. Phys. Soc. Jpn. 76, 114715 (2007). 

13 N. Hollmann, M. W. Haverkort, M. Cwik, M. Benomar, 
M. Reuther, A. Tanaka, and T. Lorenz, New J. Phys. 10, 
023018 (2008). 

14 C. F. Chang, Z. Hu, H. Wu, T. Burnus, N. Hollmann, 
M. Benomar, T. Lorenz, A. Tanaka, H. J. Lin, H. H. Hsieh, 
C. T. Chen, and L. H. Tjeng, Phys. Rev. Lett. 102, 116401 
(2009). 

15 L. M. Helme, A. T. Boothroyd, R. Coldea, D. Prab- 
hakaran, C. D. Frost, D. A. Keen, L. P. Regnault, P. G. 
Freeman, M. Enderle, and J. Kulda, Phys. Rev. B 80, 
134414 (2009). 



K. Yamada, M. Matsuda, Y. Endoh, B. Keimer, R. J. 
Birgeneau, S. Onodera, J. Mizusaki, T. Matsuura, and 
G. Shirane, Phys. Rev. B 39, 2336 (1989). 
T. G. Perring and C. D. Frost, Neutron News 15, 30 (2004). 
R. Coldea, MSLICE: a data analysis program for time-of- 
flight neutron spectrometers (2004). 

J. T. Lewandowski, R. A. Beyerlein, J. M. Longo, and 

R. A. McCauley, J. Am. Ceram. Soc. 69, 699 (1986). 

R. A. Mohan Ram, P. Ganguly, C. N. R. Rao, and J. M. 

Honig, Mat. Res. Bull. 23, 501 (1988). 

T. Kajitani, S. Hosoya, K. Hiraga, and T. Fukuda, J. 

Phys. Soc. Jpn. 59, 562 (1990). 

C. Kant, T. Rudolf, F. Schrettle, F. Mayr, J. Deisenhofer, 
P. Lunkenheimer, M. V. Eremin, and A. Loidl, Phys. Rev. 
B 78, 245103 (2008). 

G. L. Squires, Introduction to the Theory of Thermal Neu- 
tron Scattering (Dover Publications, Mineola, NY., 1996). 
A. C. Walters, T. G. Perring, J. S. Caux, A. T. Savici, 

G. D. Gu, C. C. Lee, W. Ku, and I. A. Zaliznyak, Nature 
Physics 5, 867 (2009). 

R. Coldea, S. M. Hayden, G. Aeppli, T. G. Perring, C. D. 
Frost, T. E. Mason, S. W. Cheong, and Z. Fisk, Phys. 
Rev. Lett. 86, 5377 (2001). 

N. S. Headings, S. M. Hayden, R. Coldea, and T. G. 
Perring, arXiv: 1009.2915. 

Y. J. Kim, A. Aharony, R. J. Birgeneau, F. C. Chou, 
O. Entin-Wohlman, R. W. Erwin, M. Greven, A. B. Harris, 
M. A. Kastner, I. Y. Korenblit, Y. S. Lee, and G. Shirane, 
Phys. Rev. Lett. 83, 852 (1999). 

H. M. R0nnow, D. F. McMorrow, R. Coldea, A. Harrison, 

I. D. Youngson, T. G. Perring, G. Aeppli, O. Syljuasen, 
K. Lefmann, and C. Rischel, Phys. Rev. Lett. 87, 037202 
(2001). 

N. B. Christensen, H. M. R0nnow, D. F. McMorrow, 
A. Harrison, T. G. Perring, M. Enderle, R. Coldea, L. P. 
Regnault, and G. Aeppli, Proc. Natl. Acad. Sci. U.S.A. 
104, 15264 (2007). 

T. Huberman, R. Coldea, R. A. Cowley, D. A. Tennant, 
R. L. Leheny, R. J. Christianson, and C. D. Frost, Phys. 
Rev. B 72, 014413 (2005). 



